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Atherosclerosis is a chronic inflammatory disease characterized by the accumulation of 
lipid-loaded macrophages in the arterial wall. We demonstrate that macrophage lipid body 
formation can be induced by modified lipoproteins or by inflammatory Toll-like receptor 
agonists. We used an unbiased approach to study the overlap in these pathways to identify 
regulators that control foam cell formation and atherogenesis. An analysis method inte- 
grating epigenomic and transcriptomic datasets with a transcription factor (TF) binding site 
prediction algorithm suggested that the TF ATF3 may regulate macrophage foam cell 
formation. Indeed, we found that deletion of this TF results in increased lipid body accumu- 
lation, and that ATF3 directly regulates transcription of the gene encoding cholesterol 
25-hydroxylase. We further showed that production of 25-hydroxycholesterol (25-HC) 
promotes macrophage foam cell formation. Finally, deletion of ATF3 in Apoe~l~ mice led to 
in vivo increases in foam cell formation, aortic 25-HC levels, and disease progression. These 
results define a previously unknown role for ATF3 in controlling macrophage lipid metabo- 
lism and demonstrate that ATF3 is a key intersection point for lipid metabolic and inflam- 
matory pathways in these cells. 



Cardiovascular disease is the leading cause of 
mortality and morbidity globally (World Health 
Organization Fact Sheet no. 317, http:// 
www.who.int/mediacentre/factsheets/fs317/en/ 
index.html). Atherosclerosis, a key underlying 
mechanism in cardiovascular disease, is a chronic 
inflammatory condition, characterized by focal 
thickening of artery walls and the formation of 
lipid-rich atheromatous lesions (i.e., plaques). 
These lesions contain abundant lipid-laden 
macrophage foam cells (Lundberg, 1985). The 
interaction between lipid metabolism and in- 
flainmation in plaque macrophages has been im- 
plicated in atherosclerosis progression (Moore 
and Tabas, 2011, and references therein) but 
much remains unknown about how the inter- 
action of these two processes is controlled at 
the transcriptional level. 

Many features of lipid bodies suggest that they 
represent an intersection point between macro- 
phage inflammatory and metabolic pathways. 
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Cytoplasmic lipid bodies contain a core of tri- 
glycerides (TGs) and cholesteryl esters (CEs) 
surrounded by a monolayer of phospholipids 
and proteins (Melo et al., 2011). In atheroscle- 
rotic plaque, uncontrolled uptake of oxidized 
low-density lipoprotein (oxLDL) by macrophages 
through scavenger receptors causes TG and cho- 
lesterol loading, followed by cholesterol esteri- 
fication and storage of CE in cytoplasmic lipid 
bodies (Maxfield and Tabas, 2005; Greaves and 
Gordon, 2005). Lipid body formation in macro- 
phages is also induced by a variety of inflam- 
matory stimuh, such as ToU-hke receptor (TLR) 
agonists (Nicolaou and Erridge, 2010). Although 
lipid bodies have previously been thought to sim- 
ply provide neutral lipid storage, it has recently 
become clear that they are highly regulated and 
fiinctionaUy active organelles (Imanishi et al., 2004; 
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Figure 1. Endotoxin- and lipoprotein- 
induced neutral lipid accumulations 
colocalize with lipid body marker ADRP. 

(A-C) WT BMDMs were stimulated with 
5 |jg/ml acetylated LDL (acLDL) for 4 h, stained 
for neutral lipids (BODIPY 493/503) and immuno- 
fluorescence-stained for the lipid body marker 
protein ADRP (adipose differentiation related 
protein, also known as PLIN2). ADRP detection 
was using an Alexa Fluor 633-conjugated 
secondary antibody. (For a BODIPY image of 
unstimulated macrophages, see Fig. 2 C). 

(A) BODIPY fluorescence shown in green. 

(B) Alexa Fluor 633 (ADRP) fluorescence 
shown in red. (C) Overlay. (D) Phase-contrast 
bright-field image showing refractile droplets 
that colocalize with BODIPY- and ADRP- 
positive lipid bodies. Images shown are from a 
representative of two independent experi- 
ments. Bars, 10 pm. (E) Quantitation of mean area per cell, occupied by BODIPY-positive lipid bodies in stacked confocal epifluorescence micrographs of 
WT BMDM stimulated for 18 h with 5 pg/ml acLDL, 5 pg/ml oxLDL, or 100 ng/ml LPS. Error bars represent the mean ± SEM (five fields per condition) from 
a representative experiment (two independent experiments performed). (F) Time course of mRNA levels of Atf3 in WT BMDM treated with 25 pg/ml oxLDL 
or 10 ng/ml LPS, as measured by qPCR (normalized to Eeflal]. Error bars represent the mean + SEM (two independent biological replicates). 




Bozza et al., 2007; Silva et al., 2009). The observation that key 
enzymes involved in cholesterol metabolism and fatty acid 
synthesis localize to lipid bodies suggests that they are also 
sites for hpid metabohc activity (Fujimoto et al., 2004; Zhao 
et al., 2005). Additionally, lipid bodies have been shown to 
contain eicosanoid-forming enzymes and have been demon- 
strated to be sites for eicosanoid generation (SUva et al., 2009). 

ATF3 (activating transcription factor 3), a member of the 
CREB family of basic leucine zipper transcription factors (TPs), 
has been shown to be involved in both inflammatory and meta- 
bohc pathways (Gilchrist et al., 2006; Zniuda et al., 2010). ATF3 
is transcriptionally induced in response to various conditions that 
are associated with macrophage foam cell fomiation in athero- 
sclerotic plaques including oxidative stress, hypoxia, cytokines, 
and ER stress (Hai et al, 1999; Okamoto et al, 2006; Chen 
et al., 2008; Semion et al., 2010). Notably, ATF3 is expressed at 
higher levels in human atherosclerotic vessels than in nonathero- 
sclerotic vessels, and ATF3 is strongly expressed in macrophages 
within atherosclerotic vessels (Nawa et al., 2002). Activation of 
macrophage TLRs has been shown to play a role in the progres- 
sion of atherosclerosis, and we have previously demonstrated 
that ATF3 is up-regulated early after TLR engagement and that 
it acts as a negative regulator of TLR induced cytokine produc- 
tion in macrophages (GUchiist et al., 2006). 

25-hydroxycholesterol (25-HC) is a metabolite of choles- 
terol that is produced and secreted by macrophages and has 



multiple effects on immune responses and cholesterol metab- 
olism in these cells. This oxysterol has been shown to suppress 
the production of IgA by B cells (Bauman et al., 2009), alter 
B cell migration (Pereira et al., 2010), and control the dilTer- 
entiation of monocytes into macrophages (Ecker et al., 2010). 
25-HC has long been known to suppress cholesterol biosyn- 
thesis (Kandutsch and Chen, 1975; Miller and Melnykovych, 
1984) and inhibit expression of genes involved in fatty acid 
and cholesterol synthesis (Trzaskos et al., 1989; Field et al., 
2001). 25-HC has also been shown to stimulate CE forma- 
tion (Brown et al, 1975; Miller and Melnykovych, 1984). 
The enzyme Ch25h (cholesterol 25-hydroxylase), which cata- 
lyzes the formation of this oxysterol, has been demonstrated 
to be transcriptionally up-regulated after activation of macro- 
phages by TLR agonists (Bauman et al., 2009; Diczfalusy et al., 
2009; Park and Scott, 2010). However, the transcriptional 
regulatory events that induce production of this enzyme after 
TLR activation are not fuUy understood. 

In the present study, we used an unbiased systems biology 
approach to determine that ATF3 is a key transcriptional 
regulator of the process of hpid body formation in macro- 
phages. We provide evidence that the effect of ATF3 on lipid 
body formation is mediated through direct binding to the 
promoter of the gene Ch25h, which leads to epigenetic re- 
pression of Ch25h transcription. We demonstrate that loss of 
ATF3 leads to increased macrophage foam cell formation 



Table 1 . National Center for Biotechnology Information GEO accession numbers for microarray datasets 



Strains 


Conditions 


Affymetrix GeneChip type 


GEO accession no. 


WT and Atf3-I- 


Unstimulated, oxLDL (24 h) 


Mouse Exon 1.0 ST Array 


GSE32358 


WT 


Unstimulated, LPS (4 h) 


Mouse Exon 1.0 ST Array 


GSE32359 


WT and Atf3-I- 


Unstimulated, LPS (4 h) 


Mouse Genome 430 2.0 Array (3' array) 


GSE32574 
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in vivo. Finally, we show that the severity of aortic root athero- 
sclerosis in a mouse model of diet-induced disease (the Apoe^'^ 
mouse; Smith and Breslow, 1997) is accelerated by loss of 
ATF3 and that this correlates with increased levels of 25-HC 
in aortas from Atf3^'^ animals relative to matched controls. 
Our data demonstrate a previously unknown role for ATF3 
in controlling lipid metabolism in the macrophage and show 
that this TF is a key intersection point for lipid metabolic and 
inflarmiiatory pathways in these cells. 

RESULTS 

Identification of ATF3 as a transcriptional regulator 
of lipid body formation in macrophages 

BM-derived macrophages (BMDMs) stimulated with either 
modified hpoproteins or the TLR agonist LPS accumulate 
large numbers of lipid bodies, as determined by confocal 
microscopy (Fig. 1 , A— D) and by quantitative image analysis 
(Fig. 1 E). To identify potential transcriptional regulators of 
lipid body formation, we undertook a global unbiased screen 
for TFs that are activated by both modified lipoproteins and 
LPS. Our strategy was as follows. First, we used microarray 
transcriptome profdmg data fromWT C57BL/6J BMDM to 
identify two groups of transcripts that are up-regulated by 
LPS treatment (Table SI, top 259 rows) and by oxLDL treat- 
ment (Table S2, top 35 rows), respectively. To narrow our 
search for cis-regulatory elements within the promoter re- 
gions corresponding to these transcripts, we took advantage 
of the observation that histone acetylation is associated with 
sites of TF binding to chromatin (Vettese-Dadey et al., 1996; 
Ramsey et al., 2010b) and that using histone acetylation data 
to guide a genomic sequence-based computational search for 
TF binding sites improves the accuracy of binding site predic- 
tion by '^40% (Ramsey et al., 2010b). Specifically, we per- 
formed a genome-wide chromatin immunoprecipitation 
assay with sequence-based detection (ChlP-seq) , with an anti- 
body to acetylated histones, in BMDM stimulated with 
either LPS or oxLDL (Table S3). Next, we tabulated all histone- 
acetylated regions within the promoters of the LPS- or 
oxLDL-responsive transcripts, and computationally scanned 
the regions forTF binding site sequence elements. This analy- 
sis identified eight sequence elements that were statistically 
over-represented within both the LPS- and oxLDL-associated 
promoter regions and that corresponded to TFs that were de- 
tected above background on the microarray (ATF3,TGIF1, 
API, LXR, GATAl, ZFP281, REST, and TALI). Only the 
gene Atf3 had strong differential expression under both LPS 
stimulation and oxLDL stimulation (data provided as indi- 
cated 111 Table 1). Quantitative PGR (qPGR) analysis of 
BMDM stimulated with oxLDL or LPS confirmed that Atf3 
is up-regulated by these stimuh (Fig. 1 F). GoUectively, these 
results suggested that ATF3 represents an intersection point 
for metabohc and inflammatory responses in macrophages by 
controlling lipid body formation. 

To test this hypothesis, we compared the neutral lipid 
content of Atf3~'~ and WT BMDM using flow cytometry 
for the neutral lipid fluorophore BODIPY. The difference in 
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Figure 2. Atf3~l~ macrophages have higher neutral lipid content 
than WT macrophages in the presence (or absence) of oxLDL (A) Flow 
cytometry fluorescence histograms of BODIPY-stained BMDM of the indi- 
cated genotypes, incubated for 24 h with oxLDLat the indicated concentra- 
tions. (B) Atf3-^- and WT BMDM were incubated for 24 h in medium 
containing 0, 5, or 25 ng/ml oxLDL BODIPY-stained median fluorescence 
levels are shown. (C and D) Confocal epifluorescence micrographs of BODIPY- 
stained BMDM of the indicated genotypes incubated without (C) or with 
(D) 25 (jg/ml oxLDLfor 24 h. Bars, 10 [im. Images shown are from a represen- 
tative of two independent experiments. (E) Quantitative assay for cellular CE 
and TG, normalized to total protein content, in BMDM of the indicated geno- 
types. Error bars represent mean ± SEM (three independent biological repli- 
cates). (F) Heat map showing relative transcript levels for 10 genes based on 
oligonucleotide microarray hybridization (exon array). Each column represents 
a single biological replicate for BMDM of the indicated genotype and condi- 
tion (no oxLDLor 25 ng/ml oxLDLfor 24 h), and each row represents a gene. 
Color indicates expression level relative to the gene's mean expression level in 
WT unstimulated BMDM. As shown in the heat map, microarray profiling of 
each genotype and condition was performed in three biological replicates. 



cell population median BODIPY fluorescence intensity (AT) 
between Atf3^'^ and WT BMDM was larger for cells incu- 
bated with 25 |ag/ml oxLDL (AJ = 216) or with 5 \xg/ ml oxLDL 
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Figure 3. Ch25h transcript level and 25-HC are up- 
regulated in Atf3~'~ versus WT macrophages. (A) Ch25h 
transcript levels in BMDM incubated with media alone or with 
25 pg/ml oxLDL for 24 h, measured by exon microarray. Error bars 
represent the mean ± SEM {n = 3). (B) Ch25h transcript levels in 
BMDM incubated with media alone or with 10 ng/ml LPS for 4 h, 
measured by 3' gene expression microarray. Error bars represent 
the mean ± SEM (n = 3 experiments for each bar, except for LPS 
Atf3-I-, for which n = 2). (C) Relative transcript levels of Ch25h 
in BMDM of the indicated genotypes incubated with no oxLDLor 
with 25 pg/ml oxLDL for 24 h, as measured by qPCR (normalized 
to Gapdh). Error bars represent the mean ± SEM [n = 3). (D) Lev- 
els of free 25-HC in cell extracts from BMDM of the indicated 
genotypes, measured by mass spectrometry with absolute quan- 
titation. Error bars represent the mean + SEM [n = 2). (E) Levels 
of free 25-HC in medium conditioned with BMDM of the indi- 
cated genotypes and conditions (nLDL = native LDL, 25 pg/ml) for 
24 h, measured by mass spectrometry. (F) BODIPY fluorescence 
histograms of BMDM treated with the indicated concentrations 
of 25-HC for 4 h. Data shown are from a representative of two 
independent experiments. (G) BODIPY fluorescence micrographs 
of WT BMDM grown for 24 h in medium alone, 25 ng/ml choles- 
terol, or 25 ng/ml 25-HC. Bars, 10 pm. Images shown are from a 
representative of three independent experiments. 



(A7 = 55) than for untreated BMDM (A/ = 38; Fig. 2, A and B). 
By confocal microscopy of BODIPY-stained BMDM, the in- 
creased fluorescence oi Atf3^'^ versus WT BMDM is seen to 
be the result of an increased amount of cytosohc hpid bodies 
(Fig. 2, C and D).To determine the nature of the neutral lipids 
that were increased, unstimulated BMDMs were assayed by 
lipidomic profiling using mass spectrometry. Overall, Atf3^'^ 
BMDM had 2.6-fold more CE and 1.4-fold more TG than 
WT BMDM (Fig. 2 E). 

To gain insight into the mechanism by which ATF3 controls 
macrophage Upid body formation, we performed transcriptome 
profihng oiAp^'^ andWT BMDM treated with oxLDL and 



compared them to untreated control cells. Based on the tran- 
scriptional data, the increased CE content oiAtf3^'^ versus WT 
macrophages could not be attributed to increased expression of 
transcripts encoding lipoprotein-recognizing scavenger re- 
ceptors {Msrl, Cd36, Scarbl, Cd68, and Co/eci2; list obtained 
from Greaves and Gordon, 2009), the rate-hmiting enzyme in 
the cholesterol biosynthesis pathway (Hmgcr), fatty acid syn- 
thesis enzymes (Acaca and Fasti), or sterol O-acyltransferase 1 
{Soati) in Atf3-'- versus WT BMDM (Fig. 2 F). Further- 
more, the increased CE content could not be attributed to 
decreased expression of cholesterol efflux-related transcripts 
(Fig. 2 ¥,Akai,Ahcgt and Pltp) in Atf3-'- versus WT BMDM. 
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Figure 4. ATF3 binds the promoter of 
Ch25h and histone acetylation at the 
Ch25h promoter is significantly increased in 
Atf3-I- versus WT macrophages. (A) ChlP- 
qPCR signal, relative to input chromatin, for WT 
BMDM under the indicated stimulus conditions 
and using the indicated antibodies (a-ATF3 or 
IgG from normal rabbit serum, NRS). (B) View of 
the Ch25h promoter region on mouse chromo- 
some 19, showing histone acetylation ChlP-seq 
data, CREB/ATF binding elements, and conser- 
vation data. ChlP-seq data tracks (vertical blue 
bars at 10 bp intervals) are from BMDM of the 
indicated genotypes and conditions (oxLDL, 
25 pg/ml for 24 h) immunoprecipitated with an 
antibody for acetyl-H4 or with NRS IgG as a negative control. All ChlP-seq tracks are in tags per million and are on the same vertical scale. The conservation 
data track (vertical green bars at 10 bp intervals, top row) are PHAST 30-way vertebrate conservation scores from genomic sequence, averaged over the flanking 
10 bp. Red vertical bars denote sequence matches to the CREB/ATF TP binding site motif The orange square indicates the locus that was assayed for ATF3 bind- 
ing, by ChlP-qPCR. The Ch25h gene is shown at the bottom of the diagram (thick purple bar); the direction of transcription is to the left. 
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ATF3 controls lipid body formation 
by epigenetic repression of Ch25h 

In the microarray study, the most strongly up-regulated 
gene in unstimulated Atf3~^~ versus WT BMDM (two-way 
ANOVA test;Table S4) was Ch25h, encoding cholesterol 25- 
hydroxylase (Fig. 3 A). Intriguingly, this gene is also strongly 
induced by LPS (Fig. 3 B). qPCR confirmed that loss of 
ATF3 leads to large increases in the level of Ch25h mRNA 
in unstimulated macrophages (25-fold) and in oxLDL- 
stimulated macrophages (14-fold; Fig. 3 C). The increased 
levels of Ch25h niRNA in Atf3^'' BMDM correlated with 
increased levels of the enzyme's product, 25-HC, as determined 



Figure 5. Loss of ATF3 increases neutral lipid con- 
tent and lipid body formation in macrophages from 
hypercholesterolemic mice, increases severity of ath- 
erosclerosis and increases the level of 25-HC in aortic 
tissue. Mice of the indicated genotypes were maintained 
on a Inigln-fat diet for 8 wk. (A) Flow cytometry fluores- 
cence histograms for BODIPY-stained resident peritoneal 
macrophages (B) Confocal epifluorescence micrographs of 
Nile red-stained resident peritoneal macrophages. Bars, 
10 \xm. (C and D) Left panels show aortic root lesion 
areas in male (C) and female (D) mice Long horizontal bar 
denotes mean, and the error bars denote SEM (each data 
mark denotes an individual animal). Middle and right pan- 
els show representative micrographs of aortic root sections, 
after Movat's pentachrome staining. Bars, 100 pm. 
(E) Total cholesterol and (F) TG measurements from serum. 
Error bars denote mean ± SEM (n = 6 for males, n = 4 for 
/Ipoe^" female, and /i = 5 for /Atf3"'"/lpoe"'" female). 
(G) Oil red 0- and hematoxylin-stained aortic root micro- 
graphs. Images shown are representative of two indepen- 
dent biological replicates. Bars, 100 pm. (H) Percentage of 
lesion section area staining positive for the macrophage 
marker Mac-2, by immunohistochemistn/. Error bars repre- 
sent the mean ± SEM (numbers of animals same as for 
cholesterol). (I) Ratio of 25-HC to cholesterol level, mea- 
sured by mass spectrometry with absolute quantitation, 
in sterols isolated from homogenized aortas (P < 0.01, 3.3- 
fold increase between /\poe^" and AtfS^^^Apoe^^^]. The 
long horizontal bar denotes the mean, and the short bars 
denote the SEM (each data mark denotes an individual 
aorta; n = 2 for Atf3^I^Apoe^l^, one female and one male; 
n = 4for Apoe^l^, two females and two males). Data 
shown are representative of two independent experiments. 



by mass spectrometry in both BMDM cell ex- 
tracts (Fig. 3 D) and in BMDM-conditioned 
medium (Fig. 3 E). These observations, together 
with previous data that 25-HC promotes cho- 
lesterol esterification in macrophages, led us to 
hypothesize that increased 25-HC levels in 
Atf3~^~ macrophages would promote foam 
cell formation. We therefore stimulated WT 
BMDM with 25 ng/ml 25-HC (which corre- 
sponds to a per-ceU amount that is approxi- 
mately equivalent to the amount by which 
25-HC levels increase in Atf3^'^ vs.WT macro- 
phages) and observed, using both flow cyto- 
metry (Fig. 3 F) and confocal microscopy (Fig. 3 G), that 
25-HC induces lipid body formation. 

To investigate whether ATF3 directly controls the level 
of Ch25h transcriptional activity in macrophages, we per- 
formed Chip using an antibody directed against ATF3. We 
found that ATF3 binds a predicted CREB/ATF binding site 
in the Cli25h promoter (Fig. 4, A and B). We have previously 
demonstrated that ATF3 binds histone deacetylase 1 (HDACl) 
in TLR-activated macrophages and that it acts as a negative 
regulator by epigenetic modification of cytokine gene promot- 
ers (Gilchrist et al., 2006). We therefore examined the histone 
acetylation within the Ch25h promoter using ChlP-seq and 
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found that, as with cytokines, the loss of ATF3 in macrophages 
increased the level of expression of this gene through in- 
creased histone acetylation at this locus (Fig. 4 B). 

Loss of ATF3 promotes atherosclerosis in ApoE l<nocl<out mice 

These findings led us to speculate that targeted disruption of 
ATF3 may lead to increased macrophage foam cell formation 
and increased atherosclerosis in mice. To test this hypothesis, 
Atf3^'^ mice were bred with the atherosclerosis-prone Apoe^'^ 
(apoHpoprotein E deficient) mouse strain, and the mice were 
maintained on high-fat diet for 8 wk, a time point at which 
aortic root lesions in this model would be expected to be 
macrophage-rich (Nakashima et al., 1994). Consistent with the 
in vitro findings described in the previous section, we observed 
greater neutral lipid content (Fig. 5 A) and increased Upid body 
formation (Fig. 5 B) in resident peritoneal macrophages isolated 
from the Atf3^'~Apoe~'~ animals versus from Apoe^'^ animals. 

Analysis of the extent of atherosclerotic plaque in the aortic 
root of the Atf3^'^Apoe^'^ and Apoc^'^ animals demonstrated 
that loss of ATF3 leads to increased severity of atherosclerosis 
after 8 wk on high-fat diet (P < 0.001; Fig. 5, C and D). There 
was no significant difference in overall serum cholesterol 
(Fig. 5 E) or in the serum TG levels (Fig. 5 F) between the 
genotypes. Aortic root micrographs revealed that neutral Hpid 
accumulation is predominantly cellular in both genotypes 
(Fig. 5 G). There was no significant difference in the macro- 
phage (Fig. 5 H), smooth muscle cell, or T cell (not depicted) 
content of aortic root sections between the genotypes. Fur- 
thermore, there were no significant differences between the 
genotypes, in the serum levels of the inflammatory cytokines 
TNF, IL-1(3, IL-6, MCP-1, MlPla, orMIPl(3 as detected by 
ELISA (unpublished data) . 

Based on our in vitro studies, we conjectured that enhanced 
local production of 25-HC by plaque macrophages might 
contribute to the atherosclerosis phenotype in Atf3^'^Apoe^'^ 
versus Apoe^'^ animals. Because 25-HC is thought to be pri- 
marily produced enzymatically in rodents (Johnson et al., 
1994), we hypothesized that elevated expression of Ch25h in 
plaque macrophages would lead to an increased ratio of 25-HC 
to cholesterol in the aortas of atherosclerotic mice. To test 
this hypothesis, hpid extracts from aortas were analyzed by 
quantitative mass spectrometry for 25-HC and for cholesterol. 
This analysis demonstrated an increase in the 25-HC/choles- 
terol ratio in the aortas of Atf3^'^Apoe^'^ versus Apoe^'^ 
animals (Fig. 5 I). 

DISCUSSION 

Atherosclerosis is a complex disease involving the dysregulation 
and interaction of inflanmiatory and metabolic pathways. The 
accumulation of lipid bodies in macrophages in the arterial 
intima, resulting in foam cell formation, is thought to be a 
key driver of atherogenesis. Although much is known about 
the mechanisms of foam cell fomiation, we lack complete 
knowledge of the TFs and Hpid mediators that modulate this 
process and how dysregulation of the dynamic interactions 
between them leads to disease. The inflammatory and metabolic 
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pathways in macrophages are regulated at a variety of levels that 
include cross talk and feedback between signaling pathways 
(Castrillo and Tontonoz, 2004). Systems biology provides a 
framework in which this complexity can be addressed (Zak and 
Aderem, 2009; Ramsey et al., 2010a). Systems approaches 
combine prior biological knowledge with global measurement 
technologies and computational methods to both reveal regula- 
tory interactions and place them in a network context (Gilchrist 
et al., 2006; Ramsey et al, 2008; Litvak et al., 2009; Zak et al., 
2011). The present work has used a novel systems approach to 
uncover candidate transcriptional regulators of the inflamma- 
tory and metaboUc pathways leading to macrophage foam cell 
fonnation. We used transcriptome profiUng to identify tran- 
scripts that are up-regulated by lipid body-inducing stimuli, and 
then identified histone-acetylated regions within the promoters 
of these transcripts that may indicate chromatin accessibility to 
regulatory factors. These regions were computationally scanned 
to identify a hst of candidate TFs, which was further filtered 
using transcriptome measurements. 

This approach implicated the regulator ATF3, which we 
have previously demonstrated to be a negative regulator of the 
macrophage transcriptional response to inflammatory stimuli 
(Gilchrist et al., 2006). The work presented here demonstrates 
that ATF3 also regulates lipid metabolic pathways in the macro- 
phage, at both basal and high levels of cellular hpid loading, 
and that it appears to function through a regulatory circuit 
involving epigenetic repression at the promoter of the choles- 
terol 25-hydroxylase gene Ch25h. Under loss of ATF3, the 
transcript level of Ch25h is highly up-regulated in macro- 
phages, in both the unstimulated state and in the presence of 
LPS or oxLDL. We note that Ch25h has also been reported to 
be transcriptionally up-regulated in Atf3~'~ versus WT mouse 
embryonic fibroblasts (Kim et al., 2010). Our findings estab- 
hsh that, in the absence of ATF3, the normal feedback mecha- 
nism that controls the Ch25h enzyme and its product 25-HC 
is lost, leading to high levels of macrophage production of 
25-HC and the accumulation of 25-HC in the aorta. On the 
phenotypic level, in the Apoe^'^ mouse model of diet-induced 
atherosclerosis, loss of ATF3 was found to accelerate aortic 
root atherosclerosis. As shown in Figs. 2 F and 5, we were 
unable to detect differences in serum Upid levels, serum cyto- 
kine levels, macrophage scavenger receptor transcript expres- 
sion, cholesterol efflux-related transcript expression, or plaque 
cellular composition (macrophages, T cells, or SMCs) that could 
account for this phenotype. Notably, although female mice 
had larger basal lesion area (consistent with previous observa- 
tions; Cahgiuri et al., 1999), loss of ATF3 led to a greater in- 
crease in lesion area in male mice than in female mice, an 
interesting area for possible future investigation. 

Our findings on the role of 25-HC are both consistent with 
and extend prior studies. 25-HC suppresses cholesterol synthe- 
sis largely through inhibition of 3-hydroxy-3-methylglutaryl- 
CoA-reductatase (Kandutsch and Chen, 1975), both through 
transcriptional repression of the Hmgcr gene and accelerated 
degradation of the enzyme (Trzaskos et al., 1989; Taylor, 1992). 
Correspondingly, our array data demonstrate that the transcript 
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level for Hmgcr is suppressed in cells from Atf3~^~ niice (Fig. 2 F). 
Treatment of macrophages or fibroblasts with 25-HC has also 
been reported to increase the amount of cellular CE (Brown 
et al., 1975; Miller and Melnykovych, 1984), and this increase 
has been attributed to increased availability of cholesterol 
to ACAT, and not to increased ACAT activity (Miller and 
Melnykovych, 1984; Du et al., 2004). Correspondingly, we 
observed that lipid body formation and neutral lipids are in- 
creased in response to excess levels of 25-HC in. Atp~^~ macro- 
phages, and that cellular CEs in particular are significantly 
increased in Atf3^^~ versus WT macrophages (Fig. 2 E). 25-HC 
is known to be a ligand for the oxysterol-activated TF LXR 
(Shibata and Glass, 2010) and, in that context, would be 
assumed to activate transcription of genes involved in reverse 
cholesterol transport. Our results suggest that this is not its 
dominant role, as we do not observe increased transcript levels 
of known LXR targets, such as Abcal or Ahcgl, in Atf3~^~ 
versus WT macrophages (Fig. 2 F). Intriguingly, it has been 
reported that genetic knockout of the 25- and 27-HC— 
metabolizing hepatic enzyme Cyplhi (which leads to high 
serum levels of 25-HC [Bauman et al., 2009] and 27-HC [Li- 
Hawkins et al., 2000]), in the context of the Apoe~^~ mouse 
atherosclerosis model, leads to increased aortic root lesion 
area (Umetani M., D.J. Mangelsdorf, and P.W. Shaul. 2010. 
Mechanisms of atherosclerosis: new lessons firom mouse models. 
Abstr. 18427.). 

Our data implicate ATF3 as a transcriptional regulator that 
controls foam cell formation; however, given the complexity 
of the process it is likely that other factors are involved in 
regulating foam cell formation. Furthennore, it is possible that 
loss of ATF3 contributes to lesion progression through addi- 
tional pathways, such as through increased apoptosis of vas- 
cular endothehal cells (Nawa et al., 2002). One of the strengths 
of the systems biology approach is that it allows the integra- 
tion of multiple interacting networks. The current work lays 
the foundation for the expansion of the ATF3 network that, 
in turn, will enable a more comprehensive understanding of 
the regulatory circuits that control the macrophage differen- 
tiation during foam cell formation leading to the progression 
of atherosclerosis. 

MATERIALS AND METHODS 
Mice 

C57BL/6J (WT) and Apot'-^- (B6.129P2-y4poe""'f^'VJ) mice were obtained 
from the Jackson Laboratory. Atf3~^~ mice congenic to the C57BL/6 back- 
ground (backcrossed > 10 generations) were a gift from T. Hai (Ohio State 
University, Columbus, OH; Hartman et al., 2004). Atf3~^~Apoe~^~ mice 
were generated by breeding Atf3~'~ and Apoe~^~ mice. For experiments 
involving WT and Atf3~^~ macrophages, mice were weaned onto normal 
chow and used at 8—12 wk of age. For experiments involving Apoe~'~ and 
Atf3~'~ Apoe~^~ macrophages or mice, mice were weaned onto a high-fat 
diet (protein, 17.3% wt/wt; carbohydrate, 48.5% wt/wt; fat, 21.2% wt/wt; 
cholesterol, 0.2% wt/wt; Teklad; Harlan) for 8 wk. AH experimental ani- 
mals used for atherosclerosis studies were bread and maintained at the ani- 
mal facihty of the Institute for Systems Biology. All animals were housed 
and handled in accordance with a protocol approved by the Institute for 
Systems Biology's Institutional Animal Care and Use Committee. Mice were 
sacrificed by CO2 asphyxiation and tissues/ cells were obtained as described 
in the next section. Whole blood was extracted postmortem by cardiac 



puncture and serum was obtained using a serum separator vacutainer (BD). 
Serum samples were assayed for TG and total cholesterol using fast protein 
liquid chromatography. 

Cell culture 

Cells were cultured at 37°C and in 5% CO2 in a humidified mcubator. Primary 
mouse BMDMs were prepared as previously described (Gilchrist et al., 2006). 
In brief, femoral BM cells were plated on non— tissue culture— treated plastic 
in BMDM medium (RPMI 1640 containing L-glutamine [Invitrogen] , sup- 
plemented with 10% [vol/vol] FBS [Invitrogen] that had been heat-inactivated 
[30 min at 56°C], 100 lU/ml peniciUin [CeUgro], 100 |^g/ml streptomycin 
[Cellgro], and 50 ng/ml rhM-CSF [Chiron]). On day 4, the medium was 
replaced, and on day 6, adherent cells were hfted in PBS + 1 niM EDTA, 
counted, pelleted, resuspended in BMDM medium, and plated as described 
in the next section. On day 7, cells were stimulated as indicated in the Results 
section, before being used for microscopy, flow cytometry, RNA extraction, 
or lipid extraction. Stimuli that were introduced into BMDM medium 
(as indicated in the text) were as follows: LPS (from Salmonella minnesota; List 
Biological Laboratories), nLDL (Intracel), acLDL (Intracel), CUSO4-OXLDL 
(Intracel), cholesterol (Sigma-Aldrich), or 25-HC (Avanti Polar Lipids). 
Intracel-sourced Upoproteins were derived from fresh human plasma contain- 
ing EDTA, filtered at 0.45 [im, and quantitated for protein content using a 
BCA assay (Thermo Fisher Scientific). Peritoneal macrophages were isolated 
as follows. Peritoneal exudate cells (PECs) were obtained by i.p. lavage with 
8 ml of ice-cold PBS. PECs were counted, pelleted, resuspended in com- 
plete RPMI (at 4x the cell densities used for resuspending day 6 BMDM), 
and plated onto covershps or into tissue culture dishes as described in the 
next section. Cells were incubated for 2 h and then nonadherent cells were 
removed by washing 3x with PBS. FBS-supplemented medium and oxLDL 
were assayed for endotoxin using a limulus amebocyte lysate gel-clot reagent 
(Charles River) with a limit of detection of 0.030 EU/nil. Based on this 
assay, FBS-supplemented medium contained between 0.12 and 0.30 EU/ml. 
oxLDL contained between 0.1 and 0.75 EU/ml, comparable to the FBS- 
containing medium. These values are two to three orders of magnitude 
below the LPS concentrations used to TLR4-activate macrophages (100—200 
EU/ml). Heat-inactivated FBS was also assayed for the levels of free and total 
cholesterol using the Amplex Red Cholesterol Assay kit (Invitrogen). From 
these measurements, the FBS levels of free cholesterol and CE were deter- 
mined to be 7.9 and 42.5 mg/dl, respectively. In a BMDM suspension at 
standard cell density in BMDM medium supplemented with 10% FBS, cells 
are therefore supplemented with CE at ^60 pg/cell. 

Microscopy 

50 [xl of cell suspension (containing 50,000 BMDM or 200,000 PECs) in me- 
dium was dispensed onto a #1.5 coverslip in a 24-well plate. Cells were al- 
lowed to adhere for 2 h. (For PECs, coverslips were then washed as described 
in the previous section). The volume of medium was brought up to 250 [jl/weU. 
After cells were incubated with appropriate stimulation as described in 
the text, medium was removed. For single-color (BODIPY or Nile red) 
microscopy, coverslips were prepared as follows: fixed in 10% formahn for 
15 min; incubated in PBS with 10 ^M BODIPY 493/503 (Invitrogen) for 
30 min in the dark or with 1 )j.g/ml Nile red (Sigma-Aldrich) for 5 min at RT 
in the dark; and washed 3X with PBS. For two-color fluorescence microscopy 
(ADRP + BODIPY), coverslips were prepared as follows: washed 2x with 
PBS; fixed in 3% parafomialdehyde for 20 min and washed 3x with PBS; 
permeabilized by incubation in reagent diluent (PBS with 0.1 mg/ml saponin 
and 0.5 mg/ml BSA) with 0.2 M glycine for 45 min; incubated with 1:50 di- 
lution of a-ADRP polyclonal antibody (Fitzgerald) in reagent diluent for 1 h; 
washed 3X with PBS; incubated in reagent diluent with 1 :250 dilution of labeled 
secondary (anti— guinea pig Alexa Fluor 633 conjugate) and 10 pM BODIPY 
for 1 h in the dark; and washed 4x with PBS. Coverslips were mounted to 
microscope sUdes using Vectashield hard-set mounting medium (Vector Lab- 
oratories). Confocal micrographs were acquired using a laser scanning confo- 
cal microscope (TCS SP2; Leica), in 512 X 512 scan mode, at room temperature, 
using a lOOX oil-immersion objective (1.4 NA; Leica). Confocal images in 
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Fig. 1 (A~D) were acquired at 3x zoom. Confocal images in Figs. 2 and 3 
were acquired at ix zoom. Fluorophores were excited using a 488-rmi laser 
line (for BODIPY and Nile red) and a 633-nm laser line (for Alexa Fluor 
633), and emission filters were configured using the default instrument set- 
tings for FITC (for BODIPY and Nile red) and Cy5 (for Alexa Fluor 633). 
Image acquisition was performed using the LCS software (Leica). Addition of 
scale bars and merging of color channels was performed using ImageJ software 
(National Institutes of Health). 

Lipid body quantitation 

BMDMs were plated on coversHps as described in Cell culture, and incubated 
for 18 h in media containing modified Hpoproteins, LPS, or media alone, as 
described in the Results. Cells were frxed and stained with BODIPY as de- 
scribed in Cell culture. For each coverslip, five fields were imaged using the 
confocal microscope as described in Cell culture (with ix zoom, -"^7.5 cells/ 
field), except that BODIPY fluorescence and bright-field images were ob- 
tained in stacks of 20 Z-positions over the range of positions at which the cells 
were in focus (5—12 \xm). Cells were segmented using the variance of the 
bright-field image through the Z-stack, as described previously (Selinunmii 
et al., 2009). Cellular lipid bodies were detected from the stacked fluorescence 
images, ehminating lipid bodies which are redundant in the stack, using a 
custom algorithm implemented in the CeUProfiler software (Carpenter et al., 

2006) , as described previously (Ruusuvuori et al., 2010). The total pixel area 
of lipid bodies within a field was computed, and then divided by the number 
of cells to obtain the mean hpid body pixel area per cell. 

Flow cytometry 

500 [i\ of a suspension of macrophages (3 x 10^ BMDM or 1.2 x 10^^ PECs) 
in medium was plated per well into a 12-well, iion-TC— treated tissue culture 
plate. For PECs, cells were allowed to adhere for 2 h and were washed as 
described in Cell culture. Cells were stimulated as described in the Results 
section. Medium was removed and cells were lifted by incubation in PBS + 
1 mM EDTA. Cells were pelleted, resuspended in 50 f.il of 10% formalin, 
incubated for 15 min at 4°C, pelleted, resuspended in FACS buffer (PBS 
with 2% FBS and 2 mM EDTA) containing 10 ^lM BODIPY, and incubated 
for 30 min in the dark. Cells were washed in FACS buffer, pelleted, and re- 
suspended in FACS buffer (PBS with 2% FBS and 2 iiiM EDTA) with 2.5% 
formahii. Up to 50,000 events were detected using the FACSCahbur (BD) 
with excitation using the 488-nni laser and acquisition enabled for forward 
Hght scatter, side scatter, and fluorescence emission (530/30 nni filter). The 
CellQuest Pro software was used for data acquisition. Data analysis was per- 
formed using Flowjo software (version 8.8.7; Tree Star). Events were gated 
on forward and side scatter to select the main population of cells, and the 
counts distribution for the 530/30 fluorescence emission charmel was gener- 
ated for the gated events. 

25-HC analysis 

For analysis of 25-HC in aortas, whole aortas were dissected immediately 
postmortem and homogenized in 1 ml PBS in a TissueLyser II (QIAGEN) 
for 6 mill at a frequency of 30 Hz. The homogenate was transferred to a new 
tube, the honiogenization tube was washed with 600 \x\ PBS, and the com- 
bined volume was used as the starting point for lipid extraction. For cells and 
cell-conditioned media, on day 6, 22.5 ml of BMDM medium containing a 
suspension of 16 million BMDM was plated into a 15-cm Petri dish. For 
analysis of 25-HC in cells, on day 7, cells were lifted, pelleted, and resus- 
pended in 1.6 ml PBS before addition of chloroform/methanol. For analysis 
of 25-HC in cell-conditioned media, nLDL was added on day 6, cells were 
incubated for 24 h, and 1.6 ml of medium was sampled for lipid extraction. 
For each sample type (aorta, cells, or medium), 1.6 ml of aqueous sample was 
added to 6 ml 1:2 (vol/vol) chloroform/methanol spiked with 500 pmol d3 
25-HC (Avanti Polar Lipids, Inc.) and (in the case of aortas) 5 nmol of day 
7 cholesterol (Avanti Polar Lipids, Inc.). Bulk lipid extraction and sterol iso- 
lation were performed according to a published protocol (McDonald et al., 

2007) . Sterol samples were evaporated under Nt gas using the N-EVAP ap- 
paratus (Organomation), dissolved in 50 [iX of ethanol, and then fractionated 



using reverse-phase HPLC on an Agilent 1100 system according to a pubhshed 
protocol (Rodriguez, 2004), which was modified as follows: an XTerra RP-18 
column (2.1 X 150 mm, 3.5 p-m particle size; Waters) was used, with an 
injection volume of 50 [x\ and a flow rate of 0.25 ml/min. Acetonitrile 
(ACN) was used as the organic solvent and the aqueous buffer was H2O 
containing 1 mM H3PO4. The gradient schedule was as follows: 0—5 min, 
65% ACN; 5-25 min, ramp from 65% to 100% ACN; 25-35 min, 100% 
ACN; 35-45 mm, 65% ACN. The retenUon times of 25-HC and choles- 
terol were determined to be '^ll and 25 min, respectively. For each analyte 
(25-HC and cholesterol), two fractions were collected in two 2-min bins 
(each) bracketing the analyte's retention time. The SpeedVac-evaporated 
sterol fractions were redissolved in chloroform/methanol and infused at 
2 f.il/min into a QSTAR (Applied Biosystems) or a QTOF (Waters) config- 
ured for electrospray ionization mass spectrometry (ESI-MS) in positive ion 
mode. Unmodified 25-HC was detected with (parent, fragment) ion m/z 
values (420.4, 367.4); cholesterol was detected at m/z values (404.4, 369.4). 
For each analyte, a five-point calibration curve was obtained from a serial 
dilution of standard, with a hnear fit to the caHbration data on log-log scale. 
To correct for differences in extraction efficiency between analytes, the 
abundance level for each analyte was normaUzed to the abundance level for 
the isotope-labeled standard within each biological sample. 

Lipidomic analysis of CE and TG 

BMDM of the indicated genotypes were incubated for 7 d and lifted as 
described in Cell culture. Cells were pelleted and resuspended in PBS, and 
then lysed using ceramic bead disruption. Cell lysate was assayed for total 
protein content using the BCA method (Thermo Fisher Scientific). A volume 
containing 600 pmol of 19:0 cholesterol ester (Lipid MAPS ID LMST01020002; 
Avanti Polar Lipids, Inc.) and 30 pmol of each of eight day-5— labeled TAG 
internal standards (LM-6000; Avanti Polar Lipids, Inc.) was added to a volume 
of lysate corresponding to a fixed amount of total protein (123 |ag). Total 
lipid was extracted using a modified Folch method (Folch et al., 1957). Neutral 
lipids were obtained from the total lipid extract by solid-phase extraction 
according to a previously published method (Krank et al., 2007). In brief, the 
total lipid extract was evaporated under N2 gas, redissolved in 1 ml iso- 
octane/ethyl acetate (75:25 vol/vol), sonicated, and then loaded onto sihca car- 
tridges (Isolute sflica cartridges, size XL; Biotage). Neutral lipids were eluted 
with isooctane/ethyl acetate (75:25 vol/vol), evaporated under N2 gas, and 
dissolved in 1 ml chloroform/methanol (1:2 vol/vol) containing 7.5 mM 
ammonium acetate. Neutral hpid samples were assayed for CE and TG spe- 
cies using ESI-MS in precursor-ion and neutral-loss scanning modes on a 
4000 QTRAP (Applied Biosystems) hybrid quadropole/Hnear ion trap mass 
spectrometer (Apphed Biosystems), at an infusion rate of 5 |jl/min with a 
40-min run time. Data were analyzed for total CE and TG signal using 
LipidView software version 1.1 (AB Sciex). 

Histopathology 

Aortas were flushed with PBS, fixed in 10% formaUn, embedded in paraffin, 
and the aortic root was sectioned at 4 [mi thickness. Every 10th section was 
stained with Movat's pentachrome for lesion area quantitation. Additionally, 
for each animal, three serial sections were stained for Mac-2 (macrophage 
marker), SM a actin (smooth muscle cell marker), and CD3 (T cell marker), 
respectively. Iiimiunostained sections were then stained with HRP-conjugated 
secondary and counterstained with hematoxylin. Images were acquired 
on a microscope (BX 50; Olympus) with a lOX Plan objective (Olympus) 
and a camera (EOS 5D Mark II; Canon). Quantitation was performed in 
Image Pro Plus software (6.0; Media Cybernetics). Scale bars were added 
using ImageJ. Lesion area was determined as a mean of the five sections per 
aortic root. For visualizing neutral lipids in tissue sections, forniahn- fixed, 
10% sucrose-infused, OCT-embedded frozen sections were stained with Oil 
red O and counterstained with hematoxylin. Images were obtained as 
described for the immunostained tissue sections. 

Microarray analysis 

Affymetrix exon array. Day 6 BMDMs from female mice of the indicated 
genotypes were plated at 10^ cells/well into 6-well tissue culture-treated 
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plates, in BMDM medium. On day 7, cells were stimulated as indicated, and 
total RNA was isolated using TRIzol (Invitiogen). RNA quality was ana- 
lyzed with a 2100 Bioanalyzer (Agilent Technologies). Sample niRNA was 
amphfied and labeled with One-Cycle Eukaryotic Target Labeling Assay pro- 
tocol and reagents (Affymetrix) . Labeled cDNA was hybridized to a GeneChip 
Mouse Exon LO ST array (as indicated in the Results section; Affymetrix) 
using standard protocols and reagents (Affymetrix). Imaging was performed 
with the GeneChip Scanner 3000 (Affymetrix) and the images were pro- 
cessed into probe-level intensities using GeneChip Operating Software 
(Affymetrix). Probe intensities were quantile-normalized and probe set— 
summarized (separately using exon-level and transcript-level probe sets) using 
the PLIER algorithm (Affymetrix) as implemented in the Power Tools soft- 
ware suite (Affymetrix). The probe set summarization was performed using 
probe— to— probe set mappings from the CustomCDF Project based on the 
Ensembl Exon and Ensembl Transcript databases (CustomCDF Version 13). 
Exon-level present/absent calls were made using the detection above back- 
ground (DABG; Affymetrix) algorithm with ap-value cutoff" of 0.01. Transcript 
present/absent calls were made by requiring that at least 50% of the exons for 
that transcript have a present call (See the Affymetrix technical note Identify- 
ing and vahdating alternative spUcing events: an introduction to managing 
data provided by GeneChip® exon arrays; http://media.affymetrix.com/ 
support/technical/technotes/id_altsphcingevents_technote.pdf}. For the heat 
map in Fig. 2 F, for each gene, a representative transcript was selected that had 
the highest mean log2 transcript-level probe set intensity within any of the 
four conditions. AH microarray datasets are available through National Center 
for Biotechnology Information GEO, as indicated in Table 1. Accession nos. 
for the datasets are GSE32358, GSE32359, and GSE32574. 

Affymetrix 3' array. The 3' microarray data files used to generate Fig. 3 B 
were a subset of the microarray data files used in (Gilchrist et al., 2006) cor- 
responding to WT and Atf3~^~ unstimulated and 4-h LPS-stimulated BMDM. 
The probe-level data files were reanalyzed using up-to-date genome annota- 
tions. Probe intensities were background-adjusted, quantile-normalized, and 
probe set— sunmiarized according to the Robust Multichip Average procedure 
(Irizarry et al., 2003), using the affy package (version 1.22.1) in the Biocon- 
ductor software package, and using the Ensembl Gene— based probe set mapping 
file (version 13, based on Ensembl Mm_60) from the CustomCDF project 
(Dai et al.,2005). 

Immunoprecipitation and ChlP-seq 

BMDMs were cultured from female mice of the indicated genotypes and, on 
day 7, cells were stimulated with oxLDL as described in the Results section. 
Immunoprecipitation (IP) was performed as described in Gilchrist et al. (2006), 
except that Protein A Dynabeads (Invitrogen) were used in place of agarose 
beads. Antibodies used for IP were as follows: acetyl-H4 IP (rabbit polyclonal 
IgG; Milhpore); and ATF3 (rabbit polyclonal IgG; Santa Cruz Biotechnology, 
Inc.). For ChlP-qPCR at the Ch25h promoter, the primers used were as 
follows: forward primer, 5'-TAGCAGCCCATGCTGAGACTATGT-3'; 
reverse primer, 5'-TTCTTTAGCAGGGAAAGGGAGGTG-3'; and probe, 
56-FAM/AGTGACGTCACCCTTCCGTCTCCTT/36-TAMSp (IDT). 
From the acetyl-H4 IP, a sequencing library was derived using the Illumina 
reagent kit, and 36-cycle, single-ended sequencing was performed on a 
Genome Analyzer II (Illumina). ChlP-seq data processing, mapping of reads 
to the mouse genome assembly (NCBI37), and peak detection were performed 
as described previously (Ramsey et al., 2010b). 

Combined analysis of transcriptional and ChlP-seq data 

Among Ensembl transcripts called as present in WT BMDM, differential 
expression testing was pertomied between the sample groups LPS versus no 
stimulation and between the sample groups oxLDL versus no stimulation. 
Any transcript for which the p-value for differential expression (see Statistical 
analysis) was >0.01 was excluded fr-om further analysis. For the combined 
analysis with ChlP-seq data, any transcript for which the fold change was <8 
(for LPS) or <2 (for oxLDL) was excluded from further analysis. An additional 
set of 500 background transcripts was randomly selected from a list of transcripts 



that are expressed in BMDM in any of the following three conditions: un- 
stimulated, LPS treated, or oxLDL treated. Promoter regions within ±1 kbp 
of the transcription start site of these LPS— or oxLDL— up-regulated transcripts, 
along with the background transcripts, were obtained from Ensembl (Mm_60). 
All histone-acetylated (HAc) locations within these promoter regions were 
identified by taking the base pair— level intersection of the promoter regions 
with the set of all regions identified as histone-acetylated in WT macrophages 
in any of the conditions (unstimulated, 4-h LPS-treated, or 24-h oxLDL- 
treated), using the University of California, Santa Cruz (UCSC) genome 
browser (NCBI37 assembly). These regions were extended symmetrically 
so that each region was at least 400 bp. Using the UCSC genome browser, 
the resulting three sets of mouse genomic regions (HAc regions within the 
promoters of background genes, HAc regions within the promoters of 
LPS— up-regulated genes, and HAc regions within the promoters of oxLDL— 
up-regulated genes) were mapped to sequence files for TF binding site 
motif analysis. 

TF binding site motif analysis 

Sequence files were scanned using the position-weight matrix scanning pro- 
gram Clover (Frith et al., 2004), and using all vertebrate TF binding site ma- 
trices from the TRANSFAC Professional 2010 database (BioBase). Separate 
statistical tests were performed to identify TF binding site motifs that were 
statistically overrepresented (P < 0.001, based on the empirical p-value esti- 
mate from the Clover program) within the LPS-derived sequences (vs. within 
background sequences) and within the oxLDL-derived sequences (vs. back- 
ground). TF binding site motifs that were overrepresented within both the 
LPS- and oxLDL-derived sequences were tabulated, corresponding to 10 unique 
motifs (Fig. 1 F). 

cDNA and qPCR 

Total RNA was isolated with TRIzol reagent (Invitrogen), DNase-treated 
using the Turbo DNase kit (Ambion), reverse transcribed using Superscript 
II (Invitrogen), and analyzed by real-time PGR with TaqMan Gene Expres- 
sion assays (Applied Biosystems) in a 10-^.1 reaction volume, on a 7900 Fast 
Real-Time PGR system (Applied Biosystems). TaqMan assays that were 
used are: AtfJ (Mm00476032_nil), at25h (Mni00515486_sl), and Gapdh 
(Mm99999915_gl). For Eeflai, the primers used were as follows: forward 
primer, 5'-GCAAAAACGACCCACCAATG-3'; reverse primer, 5'-GGC- 
CTGGATGGTTCAGGATA-3'; and probe, 56-FAM/CACCTGAGCA- 
GTGAAGCCAG/36-TAMSp (IDT). Data were nomialized to the expression 
levels Eeflai or Gapdh niRNA (as indicated in the text) within individ- 
ual samples. 

Statistical analysis 

Differential expression testing for microarray responses of WT BMDM to 
LPS or to oxLDL was performed using a two-tailed Student's t test on the 
log2 probe set intensities for hybridizations from each sample group, using 
Excel (Microsoft). Differential expression testing for effect of ATF3 KO on 
responses of BMDM in the presence or absence of oxLDL was performed on 
the log2 intensities using a two-way ANOVA test in R, and the p-value cut- 
oft" at an estimated gene-level false discovery rate of 0.1 (P < 2.26 x 10~^) 
was determined using the Benjamini-Hochberg procedure (Benjamini and 
Hochberg, 1995). Aortic root lesion areas were analyzed using a two-way 
ANOVA test with sex and genotype as factors, using Prism 5 software 
(GraphPad Software). 

Online supplemental material 

Lists of transcripts that are dift~erentially expressed (at least twofold) under 
LPS stimulation or under oxLDL stimulation in WT BMDM are provided 
as Tables SI (LPS) and S2 (oxLDL). ChlP-seq peak locations for acetyl-H4 
(in unstimulated, LPS-stimulated, or oxLDL-stimulated BMDM) are pro- 
vided in a UCSC BED format file as Table S3. A list of transcripts that are 
differentially expressed (at least twofold) in Atf3~^~ versus WT BMDM 
(FDR < 0.1) is provided in Table S4. OrJine supplemental material is avail- 
able at http://vAvw.jem.org/cgi/content/full/jeni.20111202/DCl. 
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